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Abstract 

A new approach is described for evaluating fracture in composite structures. This approach is 
independent of classical fracture mechanics parameters like fracture toughness. It relies on computational 
simulation and is programmed in a stand-alone integrated computer code. It is multiscale, multifunctional 
because it includes composite mechanics for the composite behavior and finite element analysis for 
predicting the structural response. It contains seven modules; layered composite mechanics (micro, 
macro, laminate), finite element, updating scheme, local fracture, global fracture, stress based failure 
modes, and fracture progression. The computer code is called CODSTRAN (Composite Durability 
Structural ANalysis). It is used in the present paper to evaluate the global fracture of four composite shell 
problems and one composite built-up structure. Results show that the composite shells. Global fracture is 
enhanced when internal pressure is combined with shear loads. The old reference denotes that nothing this 
comprehensive has been added since then. 

1.0 Introduction 

The global fracture behavior of fiber composite structures has become of increasing interest in recent 
years, because of the multitude of benefits that composites offer in practical engineering applications such 
as lightweight airframes, engine structures, space structures, marine and other transportation structures, 
high-precision machinery, and structural members in robotic manipulators. Composite structures lend 
themselves to tailoring to achieve desirable characteristics such as a high strength to weight ratio, 
dimensional stability under extreme thermal and hygral fluctuations, and the capability to allow controlled 
detectability such as in the Stealth technology. Because of the numerous possibilities with material 
combinations, composite geometry, ply orientations, and loading conditions, it is essential to have a reliable 
computational capability to predict the behavior of composites under any loading, geometry, composite 
material combination, and boundary conditions. A computational capability is also essential to design 
effective experiments for the further development of composite micromechanics theories, and to utilize 
existing experimental results in the most productive manner. In summary, the development of reliable 
computational simulation methods is necessary for the commercial maturation of composites technology. 

The behavior of composites during progressive fracture has been investigated both experimentally and 
by computational simulation (Ref. 1). Recent additions to the computational simulation have enabled 
monitoring the variations in structural properties such as natural frequencies, vibration mode shapes, and 
buckling modes during progressive fracture (Ref. 2). Existing computational capabilities in the simulation 
of structural damage and fracture of composite structures have been implemented in the CODSTRAN 
(Composite Durability STRuctural Analysis) computer program (Ref. 3). The ICAN (Integrated 
Composite ANalyzer) and MHOST computer codes (Refs. 4 to 6) are coupled to form CODSTRAN. The 
description herein is mainly to show what can be done by progressive structural fracture. Details cannot 
be included because of space limitations by conference proceedings. However, references are cited for the 
interested readers. 
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2.0 Fundamental Concept 

It is instructive to briefly describe the fundamental concepts on the origin of CODSTRAN and the 
related concepts. The most obvious one is that classical fracture mechanics are not applicable to 
composite structural fracture. A physical consideration on how a structure will fracture is very important 
in describing any new approach. 

It starts with the physical concept that a structure does not fail instantaneously, but requires some time 
for local damage to accumulate prior to structural fracture. Then the following must be happening during 
the process of damage accumulation. The process is depicted schematically in Figure 1. With Figure 1 as 
a background, this new approach to damage tolerance, structural fracture and risk is based on the 
following concepts: 

(1) Any structure or structural component can tolerate a certain amount of damage propagation 
(damage tolerance) prior to global structural fracture. 

(2) During damage accumulation, the structure exhibits progressive degradation of structural 
integrity (damage tolerance) as measured by global structural response variables such as loss in 
frequency, loss in buckling resistance or excessive displacements. 

(3) The critical damage can be characterized as the amount of damage beyond which the structural 
integrity and damage tolerance degradation is very rapid, induced by either (1) small additional 
damage or (2) small loading increase. 

(4) Structural damage tolerance degradation is characterized by the following sequential stages: 

(1) initiation, (2) growth, (3) accumulation, (4) stable or slow propagation (up to critical stage), 
and (5) unstable or very rapid propagation (beyond the critical stage) to collapse. 

The global response variables are accurately predicted by finite element structural analysis (Ref. 6). The 
composite degradation is simulated by composite mechanics in ICAN (Ref. 4). The progressive degradation 
is predicted by an incremental approach with damage tracking and the corresponding stress redistribution. 
All these are included in CODSTRAN as is depicted schematically in Figure 2 for the CODSTRAN 
computational simulation of progressive fracture. It is interesting to note that at the bottom of Figure 2 the 
constituent composite properties are degraded based on environmental and load history effects. 



Figure 1. — Overall CODSTRAN simulation. 
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Figure 2. — CODSTRAN simulation cycle. 


2.1 Approach 

The method of solution is depicted schematically in Figure 2. It starts at the micromechanics scale as 
is illustrated in Figure 3 where the local stresses and strains are described by appropriate equations 
programmed in (Ref. 5). Then, it continues to micromechanics of coupling a single fiber with a matrix 
with appropriate equation also programmed in (Ref. 5). It progresses to a mono fiber laminate; to a 
multiply laminate; to a finite element model of the structure; to structural scale that includes the boundary 
conditions, the environmental conditions, the loading conditions, and also the structural finite element 
model (Ref. 7). The scales modeled span from constituent scale to structural scale which definitely is 
multiscale. The left side of Figure 2 is “called” the composite synthesis to the single finite element node. 
The ICAN code is run for each ply in the laminate and saved so that each finite element node has its own 
ICAN code run to expedite decomposition. These ICAN runs are saved for the downward composite 
structural decomposition as noted in the right side of Figure 2. Once the finite element solution of the first 
load increment with internal forces/displacements has been obtained at each node then the downward 
decomposition starts. It is noted that the finite element solution requires nodal information because it is 
computationally more expedient for the composite decomposition to be performed (Ref. 8). Then the 
decomposition is repeated by using the ply information stored in the synthesis process. The mono ply 
stresses/strains are then evaluated by using the schematic in Figure 3 where the local failures are 
identified. If any failures occurred at this level, the respective stiffness and fractured region are eliminated 
for the second simulation. The process continues until local incremental convergence has occurred. At 
this point the load is increased by the second increment. Loading increments are progressively larger at 
the beginning until local fracture is detected. Then the load increment is reverted back to the last 
increment and is progressively halved until convergence is achieved and the next load increment is 
applied by a value equal to the previous load increment. Figure 4 illustrates this concept. Therefore, the 
solution is incremental from the micromechanics scale to the structural local/global convergent scale. The 
structural dynamics equations solved by the finite element in CODSTRAN, which have global variable 
convergence criteria, are summarized in the chart as is depicted in Figure 5. These equations are solved at 
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The six ply stress components are along the material axes. 

S L11C < a Lll < ^HT 

S L22C < a L22 < S L22T 
S L33C < ct L33 k S L33T 
S L12(-) < a L12 < S L12(+) 
s L23(-) < a L23 < S L23(+) 

S L1 3(-) < a L13 < S L13(+) 



MDE combined stress failure criterion. 

F = l-[ (^Liia/SLija) 2 + (^L22b /S L22b> 2 ‘ K L12ab( a Llla /S Llla > < c L22b /s L22b) 
+ (a L12S /S L12S^ !• 

Where if 1 > F > 0 No ply failure 

F = 0 Ply failure imminent 

1 < F Ply failure has occurred 

Figure 6. — Ply failure criteria, maximum stress failure criterion. 


each load increment. There is another level of ply failure criteria. This is a stress failure criterion with a 
combined stress failure function as depicted in Figure 6. The combined stress failure criterion is applied 
first. Then ply dominant stress is identified and that stress is used in the constituents to identify which 
region has failed. Therefore, the solution is robust and quite reliable as will be illustrated subsequently by 
the solution results of the sample cases. 


3.0 Sample Cases Results and Discussions 

3.1 Composite Shell With Defects 

The first illustrative sample problem is a composite shell subjected to internal pressure and having a 
structural condition of a through-the-thickness longitudinal defect, as shown in Figure 7(a) (Ref. 8). The 
type of composite and the laminate configuration of the shell are also shown in the title of Figure 7. 
Additional details are described in (Ref. 8). The environmental conditions are noted in the small table 
insert to the right of the shell. The results are plotted pressure versus damage percent in part (b) top right. 
Third vibration frequency versus pressure part (c), down left; and third vibration frequency versus 
damage percent (d), down right. Each plot has six different curves, one each for environmental effects. 
The very top curve (o) is with no environmental effects. The second from the top curve (□) is room 
temperature and one-percent moisture content by volume. The third from the top curve (A) is for the 
temperature 99.9 °C (200 °F). The fourth from the top curve (0) represents the combined temperature 
moisture effects 200 °F with one-percent moisture by volume. The fifth from the top curve (n) is for 
temperature 149 °C (300 °F) only. The last curve (V) is for the combined environmental effects 148.9 °C 
(300 °F) with one-percent by volume moisture. Note that the 148.9 °C (300 °F) temperature only curve 
shows the second greatest failure pressure. The reason is that that shell has the lowest residual stress that 
counteracts the temperature degradation effects. The important point to observe in these results is that the 
environmental effects have substantial structural integrity degradation effects. The curves plotted in 
Figure 7(c) show the significant degradation on the third vibration frequency. The structural degradation 
effects are also significant when the third vibration frequency is plotted versus damage percent. 
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Figure 7. — Hygrothermal effects — T300/epoxy [90 ±1 5] s shell with fiber volume 
ratio 0.6. (a) Geometry and environment, (b) Pressure, (c) Vibration frequency, 
(d) Buckling load (1 psi = 6.9 Pa). 


3.2 Composite Shell With Defects and External Pressure 

The second illustrative example is a composite cylindrical shell under external pressure with no 
defects, and with two types of defects — surface plies, and mid-thickness (Fig. 8) (Ref 9). The shell 
composite type and laminate configuration are shown at the bottom of the figure as well as defect sizes. 
The shell with the surface defects showed limited damage tolerance of about 149 Pa (50 psi). The shell 
with the mid-thickness defects have rather insignificant effect compared to defect free shell which 
exhibited no damage tolerance at all. The important conclusion from the evaluation of this shell is that 
only surface defects show a limited damage tolerance when the shell is subjected to external pressure. The 
reason is that compression has an enhancement in small defects. Additional details for this composite 
shell are described in Reference 9. 

3.3 Composite Thick Shell 

The third illustrative example is a thick shell subjected to external pressure and with three types of 
defects as shown in Figure 9, where the composite type, laminate configuration and finite element model 
are also shown (Refs. 10 and 1 1). The location and shape of the defects are shown in Figure 10. Results 
obtained are shown in Figure 11. This figure is very important because it shows the damage initiation and 
fracture of the defect free and defected shell. It is interesting to note the ply lay-up of the shell thickness 
in Figure 10. The damage initiation is shown as ratios of the defect free shell which is unity (1.0). The 
damage initiation of the defect free shell is 0.84. This value indicates that the shell has 16 percent damage 
tolerance from initiation to global fracture. The mid-thickness defects exhibited a 0.75 initiation and 0.77 
global fracture. The damage tolerance of the shell with the mid-thickness defects had almost no damage 
tolerance. That is the shell exhibited a rather brittle behavior. The shell with the inner surface defects had 
a 0.45 damage initiation and 0.85 global fracture. This shell had the greatest damage tolerance of 
40 percent from its initial damage to its global fracture. For additional details see Reference 11. 
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Figure 8. — Thick shell fracture under external pressure. 



Composite Shell T300/Epoxy[902/±15/90 2 /±15/90 2 /Tl5/90 2 ] 
shell diameter = 40 in. length = 80 in. 

612 nodes, 576 quadrilateral elements 
initial fiber defect in 2 adjacent hoop plies 
defect extends 5 in. along axial direction of shell 


Figure 9. — Shell structure evaluated (1 in. = 1.5 cm). 
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Figure 10. — Shell laminate structure schematic indicating initial defects 
(1 in. = 1.5 cm). 



Figure 11 . — Summary of results. Composite shell 
T300/epoxy [90 2 /±1 5/90 2 /±1 5/90 2 /+1 5/90 2 ]. 


3.4 Composite Thick Shell Global Degradation 

The fourth illustrative example is another thick shell under external pressure as shown in Figure 12 
where the shell composite system, laminate configuration, finite element model are also shown. This shell 
was analyzed for degradation in the frequency and the buckling load as the damage propagated along the 
longitudinal direction as shown in Figure 13. The buckling load did not degrade until the damage length 
was about 42 cm (28 in.) long. After that the buckling load degradation was relatively great with global 
collapse of about 80 psi down from 2.3 KPa (340 psi) or a degradation of about 1.8 KPa (260 psi). This is 
a very interesting result because it indicates the buckling of a composite thick shell has a relative large 
damage tolerance with respect to buckling resistance. Though frequency degradations are not shown here, 
these degrade slower than the buckling load. Additional details are described in Reference 1 1 . 


NAS A/TM— 20 12-217680 


350 q 



T300/Epoxy Composite Structure 
Fiber volume ratio is 0.6 
[[[[[[90 2 /±15/90 2 /Tl5/90],]JJJJ, 

800 nodes, 768 quadrilateral elements 
First Natural Frequency = 25.6 Hz. 
First Buckling Load = 337 psi. 
service load = 200 psi. 


Figure 12. — Composite shell: 576 plies 
(2.88 in. thick) diameter = 15 ft; length = 50 ft 
(1 psi = 6.9 Pa; 1 ft = 18 cm). 
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Figure 13. — Damage tolerance for composite shell 
(1 psi = 6.9 Pa). 
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Figure 14. — Stiffened composite shell panel. 


3.5 Composite Built-Up Structure 

The fifth illustrative example is for the fracture of a built-up composite structure (Ref. 12). The 
schematic of this structure is depicted in Figure 14 where the composite type and a global finite element 
model are shown of a segment of the composite structure as well as the loading conditions. The composite 
laminate configuration is shown in Figure 15 where the finite element model of the internal structure is 
also shown. The tensile results are summarized in Figure 16 where the load is plotted against damage 
percent. Figure 16 illustrates the load which initiated damage and the load at which global fracture 
occurred. The points to be observed in Figure 16 are: (1) the tensile load initial damage starts later and 
caused maximum damage propagation at fracture; (2) shear with tension exhibited the least damage 
propagation; and (3) the combination of axial tension with shear and with pressure had the intermediate 
damage propagation and tensile load to fracture. An important conclusion is that the internal pressure 
stabilizes fracture and shear load. 
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Figure 15. — Laminate structure schematic. 
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Figure 16. — Axial tension load and damage progression. 


A comparable plot is shown in Figure 17 with compressive load. In this figure it is shown that the 
combined load of compression, shear and internal pressure sustained the graded amount of damage 
propagation. However, the compression only and compression with shear caused about the same amount 
of damage propagation but global fracture occurred at different loads. Compression only sustained about 
50 percent more load at global fracture than compression combined with shear. In this case the internal 
pressure stabilized the structure to sustain greater damage tolerance at fracture. A summary of all the 
loading conditions, where the fracture initiated with each loading condition, damage propagation (growth 
plus progression), and global fracture are shown in Figure 18. The information in Figure 18 is unique and 
demonstrates the utility of the CODSTRAN computer codes as well as the wealth of information obtained 


NAS A/TM— 20 12-217680 


10 




Damage, percent 

Figure 17. — Axial compressive load and damage progression. 
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from incremental up-dating solution methods. It is interesting to check whether the incremental up-dating 
can evaluate nonlinear problems. An example of this nonlinearity is illustrated in Figure 19. In this case 
the nonlinear compression load versus applied micro-strain is compared from that obtained by a 
commercial code nonlinear analysis and that predicted by CODSTRAN. It is observed that CODSTRAN 
predicted a substantially different nonlinear curve compared to that obtained by using a commercial finite 
element computer code. Additional discussion on this problem is provided in Reference 12. Another 
important conclusion is that up-dated incremental solutions treat large nonlinearities extremely well. It is 
a general conclusion of this investigation that CODSTRAN is a unique code for treating fracture 
progression and global structural fracture in composite structures with multifunctional multiscale 
problems. Therefore, CODSTRAN is a truly “Desk-Top” portable “Testing Machine (Ref. 10).” 
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Applied microstrains 


Figure 19. — Compressive load with end displacement AS-4/HMHS[[0/±45/90] s ] 6 
(1 kip = 6.9 KPa). 


4.0 Concluding Remarks 

Up-dated computational simulation is one direct approach to evaluate fracture in composite 
structures. CODSTRAN is a stand-alone general purpose integrated multiscale/multifunctional computer 
code which consists of several modules including micromechanics through structural analysis. It is 
applicable to general classes of composite structures. Composite shells fracture investigated herein 
included defect free shells and shells with defects. The simulation results presented are from the 
microscale to global structural fracture. A built-up composite structure subjected to combined loads was 
evaluated from micromechanics fracture to global fracture. Results from all of the above problems 
indicate that shear load combined with tension or compression stabilize the solution as shown by the 
greater damage sustained at global structural fracture. Embedded defects have no influence in the global 
shell fracture when the shell is subjected to internal pressure. 
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